Long-range entanglement in the Dirac vacuum 



J. Silman and B. Reznik 
1st February 2008 

School of Physics and Astronomy, Raymond and Beverly Sackler Faculty of Exact Sciences, 
Department of Physics and Astronomy, Tel-Aviv University, Tel- Aviv 69978, Israel 



Abstract 

Recently, there have been a number of works investigating the entanglement 
properties of distinct noncomplementary parts of discrete and continuous Bosonic 
systems in ground and thermal states. The Fermionic case, however, has yet to 
be expressly addressed. In this paper we investigate the entanglement between 
a pair of far-apart regions of the 3+1 dimensional massless Dirac vacuum via a 
previously introduced distillation protocol [B. Reznik, et al., Phys. Rev. A 71, 
042104 (2005)]. We show that entanglement persists over arbitrary distances, and 
that as a function of L/R, where L is the distance between the regions and R is their 
typical scale, it decays no faster than ~ exp [— (L/i?)^]. We discuss the similarities 
and differences with analogous results obtained for the massless Klein-Gordon 
vacuum. 



Entanglement in spatially extended many body systems and quantum field theories is the focus 
of increasing attention. Part of this is directed at understanding the entanglement properties of 
noncomplementary parts of a system, such as far apart regions of vacuum 0121 El Bill] and thermal 
states [Sl m, or widely separated segments in ground states of the Bose-Hubbard Hamiltonian 



IH], chains of trapped ions (HI and harmonic oscillators QHl- In this regard relativistic vacua are 
especially interesting as they provide us with an oppurtunity to study physical systems with a well 
defined notion of locality. 

In this paper we investigate the entanglement between arbitrarily distant regions of the free 
massless Dirac vacuum. For Bosonic systems, the expansion of the vacuum in terms of two-mode 
squeezed states of oscillators residing in the two complementary spacetime wedges a; > and x < 0, 
used in the derivation of the Unruh effect, explicitly shows that the vacuum is entangled jlll I12j . 
This result is a special case of a general modewise decomposition theorem pertaining to a certain 
class of Bosonic Gaussian states JHl ■ analogous theorem exists for Fermionic Gaussian 
states 51]. (Indeed, the Unruh effect holds also in the Fermionic vacuum The state of a 

pair of noncomplementary parts of a system, however, in general is mixed, so that a modewise 
decomposition is impossible |17 . Working directly with the system's degrees of freedom, especially 
when of a great or inifinite number, proves difficult then. A most effective and relatively simple way 
to tackle this problem is the use of entanglement distillation protocols. Even though such protocols 
have proved most convenient in the study of the entanglement between abitrarily distant regions 
of the Bosonic vacuum S E], the Fermionic case has thus far not been expressly addressed. 
Using a previously introduced distillation protocol OE], we explicitly show that results analogous 
to those obtained for the Bosonic vacuum are true of the Fermionic vacuum as well, namely, that 
entanglement persists between arbitrarily far-apart regions and that as a function of the ratio of 
the separation between the regions L and their typical scale R, the entanglement decays no faster 
than ~ exp —{L/Kf 

The concepts of entanglement and locality are nontrivial in Fermionic systems. We therefore 
begin by explaining them briefiy and contrast with the Bosonic case. Suppose we have a system of 



Bosonic modes. As 



— 6ij and [a^, aj] — 0, the Hilbert space is a direct product of the Hilbert 



spaces of each of the modes. Hence, it is meaningful to consider the entanglement between different 
sets of modes, with the partition unequivocally defining locality. For Fermions |aj, a|| — Sij 
and {tti, aj} ~ 0. The Hilbert space therefore lacks an analogous direct product structure. If the 
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assignation of sets of modes to different parties is to have any meaning at all (that is, if we do 
not want to give up locality) , we must restrict the set of observables in such a way that acting on 
an arbitrary composite state of any two distinct sets of modes, say A and B, with an observable 
comprised solely of modes in A, does not change the expectation values of any observable comprised 
solely of modes in B (and vice- versa), nor increase the entanglement between the sets. Now, let 
Oa and Ob be arbitrary sums of products of even number of modes in A and B, respectively, then 
Oa, Ob = 0. It is not hard to see that this characteristic of Fermionic modes, together with 
the anti-commutation algebra of the modes, restricts the set of observables to those that can be 
constructed out of products of an even number of modes JHl ■ For pure states entanglement between 
two sets of modes is then defined as usual, i.e. a pure composite state of two sets of modes, (f), is 



entangled iff there exist observables Oa and Ob such that {OaObJ^ 7^ \^)^ \^)ip' course 
this is not true of mixed states. Indeed, except in 2 x 2 and 2x3 dimensions ^21 ED], no necessary 
and sufficient criterion to estabHsh mixed state entanglement is known, regardless of the statistics. 

Moving on to relativistic quantum field theory (QFT), the requirement of Lorentz covariance and 
that the energy spectrum be bounded from below constrains the set of possible algebras of modes 
to the familiar commutation/anti-commutation relations for Bosons/Fermions. As an example 
consider the Dirac field 

{1/., (x, t) , ^] (y, t) } = 5,j5 {x~y), {V-. (f , t) , {y, t)}^0. (1) 

The subscripts denote the spinorial indices, which together with the position x label the modes. If 
in addition we want the theory to be causal, we must require that observables be respresented by 
bilinear expressions in the fields. 

It is important to note that there is a difference in what is meant by "local" in quantum infor- 
mation theory (QIT) and QFT settings. As explained above, in QIT it is the different parties that 
define locality. However, in QFT it is causality which defines locality, i.e. an operator acting at two 
or more spaceHke related coordinates is nonlocal. In this paper, locality in the QIT sense enters via 
the assignation of causally disconnected regions to different parties, leaving us with much greater 



3 



latitude in our choice of "local" operations than that afforded by the tight constraints of QFT. 

The protocol employed in (2 E] consists of the finite duration coupling of a pair of initially 
nonentangled two-level point-like detectors to the studied field, in its vacuum state, at two different 
locations. The duration of the coupling determines the size of the regions "probed" and is taken to be 
much smaller than the distance between the detectors, which therefore remain causally disconnected. 
Under these conditions, a final entangled state of the detectors means that entanglement persists 
between the regions. A similar but suitably adjusted protocol is employed here. We therefore make 
use many of the results obtained in and forego rederivation. 

The Dirac equation is given by 

d 

{a-p + (3m)ij{x,t) = i—ii^{x,t), (2) 

where the ai and (3 are any 4x4 matrices satisfying {ai, aj} = 25ijl, 0^ = 1, and {a, /?} = 0. 
In the absence of a mass term, in the Weyl representation of the ai and /3, the Dirac equation 
decouples into a pair of equations, the Weyl equations 

d d 

a ■ piprix, t) = i—ipr{x, t) , -a ■ pipiix, t) = i—iji{x, t) , (3) 

where the ai are the Pauli matrices. The two-component fields {x, t) and ipi {x, t) describe 
right and left handed particles and anti-particles, that is quanta of positive and negative helicity, 
respectively. We can therefore begin by studying a vacuum of definite handedness, say the right 
handed vacuum. In terms of a Fourier expansion 

Op and 6^^, 6^ are the creation and annihilation operators for right handed particles and 
anti-particles, respectively, satisfying |a^, = ^^?^} ~ (^tt)^ 5 (p — (7^ with all other anti- 

commutators vanishing, while Ur (p) and Vr (p) are the corresponding two-component "spinorial" 
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coefficients and are understood to be normalised to unity. 

Due to the fact that the field is Fermionic, we must couple to field bilinears to realise the 
distillation protocol (See earlier discussion.) Perhaps the most natural choice is the field's 
charge density N (-0^ {xi, t) tpr {xi, t)), where for convenience we have chosen to normal order (N). 
Setting up a pair of two-level detectors at x^and xb, in the Dirac interaction picture the coupling 
term is given by 



(5) 



=A, B 



Here {t) governs the strength and duration of the coupling, and VLi is the energy gap of detector 



, with T and T 



i. The corresponding evolution operator is U{T/2) ~ Texp —i J'^^^^^dtHc (t) 
denoting time-ordering and the duration of the interaction, respectively. As discussed above, we 
set L 3> T (L \xb — xa\), and take the initial state of the detectors to be separable. 

Once the interaction is over, in the basis {||, i]', jj,, 11}, the partial transpose of the detectors' 
reduced density matrix is given by 



/ 1- 



PT 
PAB 



EA^-WEBf 



\ {Eb I Ea) 



Q {Ea\ Eb) ^ 

llBsf -{Q\Xab) 

(XabIO) WEa^ 

\\XABf j 



(6) 



where 



\E, 



T/2 



-T/2 



dte, (t) e'"-*iV {^Pl (x,, t) t)) |0) , 



(7) 



rT/2 

\Xab) ■■= / dtdt'eA {t) en (f) eK^^^*+"«*')7V {xa, t) 0, {xa, t)) N {^l {xb, t') (xs, i')) |0) 

J-T/2 

(8) 



Using the Peres criterion ,1.9|, we find that the detectors are entangled (i.e. that the partial transpose 
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has negative eigenvalues) if 



\{0\XAB)\'-mf\\EBf>0. (9) 



Physically speaking, this translates to the requirement that the probability of exchange of a right 
handed virtual particle - anti-particle pair between the detectors be greater than the product of 
the probabilities for the on-shell emission of a right handed particle - anti-particle pair by the same 
detector. 

For temporally symmetric window functions, a somewhat lengthy calculation shows that the 
above condition takes on the explicit form (see appendix for details) 



— Sin (llIiL) Sin {1JO2L) 



—uJiCos (wi-L) Sin (UJ2L) — oJ^Sin (uJiL) Cos {UJ2L) 



Ia {^a t^2) es {^B - LOi- L02) 



— / duju^eA {^A +0;) / dioLu^eB 
^5 Jo Jo 



> 



{nB + ujf, (10) 



where Ci is the Fourier transform of e^. 

This is to be compared with the condition obtained for the massless real Klein-Gordon field [2] 

/•oo 2 /.oo poo 

/ — 5m (wL) (ri^ + o;) (f^B — cj) > / dujui \eA {^a + ^^)\ / dojuj \eB {^b + ^)\ 
Jo L Jo Jo 

(11) 

The two conditions bear similarity. The analysis performed in shows that Eq. Ijlljl can be 
satisfied if we choose i^A — ^) such that it oscillates as Sin (uiL), that is faster than any of its 
Fourier components, over a finite integration regime |2;^ll24j . Indeed, such a choice [221 can render 
the exchange probability arbitrarily larger than the product of the emission probabilities. Suppose 
in our case we take the superoscillatory transform to oscillate Hke Cos (ujL) over a suitably chosen 
integration regime. Then the specially tailored form of the superoscillatory transform guarantees 
that the first term on the LHS of Eq. ifTIijl is much greater than the RHS. But for precisely the 
same reason it is much greater than all other terms on the LHS, and Eq. ifTTill is satisfied. It follows 
that entanglement persists between arbitrarily far-apart regions of the massless Dirac vacuum of 
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Figure 1: Spacetime at the start of the interaction t — ~T/2. A and B denote the locations of the 
detectors. It is only operators within the two spheres that contribute to the distilled entanglement. 

quanta of definite handedness, and that the lower bound obtained in P] holds here as well. That 
is, in the limit L/T ^ 1 the entanglement, quantified by the negativity J\f, scales no faster than 



exp 



{L/ry . Now as in the duration T the detectors probe a spherical region of radius R = T 



(see Fig. 1), we arrive at the aforementioned lower bound ~ exp — (L/R)'^ 

Not surpirsingly, for the left handed vacuum the condition for entanglement is identical. This 
means that double the amount of entanglement can be distilled by coupling to the total charge 
density, that is the sum of the charge densities of right and left handed quanta. 

Before we conclude, there are three questions that need to be addressed. First, as previously 
mentioned, an identical bound for the entanglement obtains for the Klein-Gordon vacuum. The 
question arises as to whether this reflects some sort of universality or is just an artifact of our 
distillation protocol [2^. Naively, it might be expected that the comparatively "poorer" structure 
of the Dirac Hilbert space, resulting from the anti-commutative nature of the fleld, leads to a faster 
decay. However, the fact that our distillation protocol is perturbative means that in the Klein- 
Gordon fleld case, the Hilbert space's full structure does not come into play. This may very well be 
the reason for the identical bound, but to say more would be pure speculation. 

Second, it is natural to ask whether the correlations giving rise to this entanglement can be 
attributed to a local hidden-variable model. In the Klein-Gordon fleld case, we were able to show 
PI the detector's flnal state exhibits "hidden" nonlocal correlations [221, the sense that after local 
flltering jJH] an EPR state can be distilled. However, the same does not true here, due to the 
presence of terms in the reduced density matrix, Eq. lO, not present in the Bosonic equivalent. 
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which prevent the distillation of an EPR state. This question therefore remains open. 

The third question is how the results obtained change in the massive case. Obviously, the 
presence of a mass term adds another scale to the problem. From the point of view of our distillation 
protocol, the Dirac equation no longer decouples and it is hard to see how use can be made of a 
superoscillating function to satisfy the resulting inequalities (without which we do not know how to 
distill entanglement at arbitrarily long distances). Nonetheless, the fact that the main contribution 
to the entanglement arises from high frequencies (see pj) suggests that for a comparatively small 
mass our results should remain unchanged [29| . 



Appendix: Details of calculations 

To clarify some of the physical content behind the condition for entanglement, Eq. l[T?1|l . we outline 
here the important steps in its derivation. 

As already noted, in the absence of a mass term, the Dirac equation decouples into a pair of 
equations for quanta of a definite handedness. However, it is only in the Weyl representation of the 
Gamma matrices 

/ n 1 \ / n rr. \ 

(12) 

that these equations reduce from four component equations to two. Taking the spinors to be 
normalised to unity we then have 



1 { Px ^ iPy 

Ur [p] = Vr [p) = , , = I , Ul{p) ^Vl{p) ^ Ur {-p) ■ (13) 











( 




70 = 

























\/'^P{P~Pz 



P-Pz 



Focusing on the right handed vacuum, in terms of the Fourier expansion of the field the emission 
and exchange terms are given by 

\m' = J (P) 4 (p) Vr (q) {n,+p + q)\' , (14) 
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{0\Xab)^ - f ^j^vl (p) Ur (q) 4 (q) Vr ip) e^(P+'r>-(-~--^-)iA (17a + P + q) iB {ns-p-q), 

(15) 

where we have ah-eady carried out the temporal integration. Plugging in the expressions for the 
spinors we get 



(Ppd^q if P ■ q 



{0\Xab)^- f ^^ifl + ^^)e^(^~+«>(-~---~-)eA(17A+P + 9)6B(f7B-p-'z). (17) 
J {2ttY 2 V pq J 

In spherical coordinates the integration over angles is straightforward. 

||iJ,f = 2 r^pV|e-(a:+P + <?)r, (18) 
JO 27r 



(0 I Xab) = -2 r \Cos Hp + q)L)- ^Cos (pL) Sin (qL) - ^Sm [pL) Cos {qL) 

Jo {2tt) ^ L L 

pq 1 
+ j^Sin {pL) Sin (qL)^ eb {^a +P + q)eB {^b - P - q) , (19) 

where L \xa ~ xb\- If we now switch to the variables uj — p + q and v — p~ q, then integrating 
over V, Eq. ltTn)l quickly follows. 

We note that the Cos {uoL) term on the LHS of Eq. ltT?1|l arises from the angular integration 
over both the ^ and ^ terms resulting from the spinor products. Were the ^ term absent, we 
would not be able to distill entanglement at any distance L. It is interesting that this implies that 
we cannot reaHse our distillation protocol in the real scalar vacuum via a square coupHng. 
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